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Abstract 



The properties of the low-energy dipole response are investigated for the proton- 
i rich doubly magic nucleus ^^Ni, in a comparative study of two microscopic models: 

I fully self-consistent Relativistic Random-Phase Approximation(RRPA) based on 

the novel density-dependent meson-exchange interactions, and Continuum Random- 
Phase Approximation(CRPA) using Skyrme-type interactions with the continuum 
I properly included. Both models predict the existence of the low-energy soft mode, 

lO ■ i.e. the proton pygmy dipole resonance (PDR), for which the transition densities and 

RPA amplitudes indicate the dynamics of loosely bound protons vibrating against 
the rest of the nucleons. The CRPA analysis indicates that the escape width for the 
proton PDR is rather large, as a result of the coupling to the continuum. 

^ ■ PACS: 21.10.Gv, 21.30.Fe, 21.60.Jz, 24.30. Cz 



One of the major challenges in the region of unstable nuclei is the under- 
5^ I standing of soft modes of excitations which involve loosely bound nucleons, 

nucleon halo, or skin. In particular, in neutron-rich nuclei, nucleons from the 
neutron skin may give rise to a low-energy dipole mode, known as a pygmy 
dipole resonance (PDR) [1]. Some experimental evidence about the low-lying 
dipole excitation phenomena in neutron-rich nuclei is available from (i) elec- 
tromagnetic excitations in heavy-ion collisions [2], and (ii) nuclear resonance 
fluorescence experiments for nuclei with moderate neutron to proton number 
ratios, e.g. ^^^'^^^Sn [3], ^^SBa [4], Pb [5], and Ca [6] isotopes, and N=82 iso- 
tones [7]. On the other side, a variety of theoretical models have in recent years 
been employed in studies of the low-lying El strength: continuum quasiparti- 
cle RPA formulated in the coordinate-space Hartree-Fock-Bogoliubov frame- 
work [8], continuum RPA with Woods-Saxon potential for the ground state 



* Corresponding author Email: nils.paar@physik.tu-darmstadt.de 
^ Permanent address: Joint Institute of Nuclear Research, Dubna, Russia 



Preprint submitted to Elsevier Science 



9 February 2008 



and Landau- Migdal force in the residual interaction [9], the time- dependent 
density-matrix theory [10], Skyrme Hartree-Fock + quasiparticle RPA with 
phonon couphng [11,12], and the quasiparticle-phonon model [3,5,13]. In the 
relativistic framework, properties of low-energy excitations have been system- 
atically studied in the relativistic RPA (RRPA) [14], and the relativistic quasi- 
particle RPA (RQRPA) [15,16,17,18]. Both experimental and theoretical stud- 
ies qualitatively agree on the global properties of the isovcctor dipole response, 
i.e., as the number of neutrons increases along the isotope chain, the transi- 
tion strength distribution is characterized by the appearance of pronounced 
low-lying El strength. Within the relativistic R(Q)RPA studies, it has been 
shown in a fully self-consistent way that the low-lying pygmy state represents 
a genuine structure effect: the neutron skin oscillates against the core, exhibit- 
ing a collective nature from medium towards heavier nuclei [14,15]. A detailed 
quantitative description of the low-lying El strength in neutron-rich nuclei is 
essential for calculations of radiative neutron capture rates in the r-process, 
and elemental abundances from nucleosynthesis [19]. 

The structure of nuclei on the proton-rich side is equally important for reveal- 
ing many aspects of the underlying many-body problem and properties of the 
effective nuclear interactions. At the same time, a quantitative description of 
proton-rich nuclei is of a particular importance for describing the rapid pro- 
ton capture process of nucleosynthesis. Proton-rich nuclei are characterized 
by unique ground-state properties such as j3 decays with large Q values, and 
direct emission of charged particles. The two-proton ground state radioactiv- 
ity has recently been observed as a decay mode of ^^Fe [20,21]. At the present 
time, knowledge on dipole excitations in nuclei towards the proton drip-line is 
rather limited. Only very recently, the first microscopic theoretical prediction 
of the proton pygmy dipole resonance has been indicated in the framework of 
RQRPA based on the Relativistic Hartree-Bogoliubov model (RHB) [22]. For 
nuclei close to the proton drip-line, the model calculations predicted the oc- 
currence of pronounced dipole peaks below 10 MeV in excitation energy, due 
to collective vibrations of loosely-bound protons against the proton-neutron 
core. 

In this letter, we employ two theoretical approaches in a comparative study 
of the low-lying El strength of the proton-rich nucleus ^^Ni: the fully self- 
consistent Relativistic RPA (RRPA) based on novel density-dependent inter- 
actions, and the Continuum RPA (CRPA) with Skyrme-type interactions. An 
essential objective of this study is to ensure that the proton PDR is inherent 
for different models, based on different assumptions and effective interactions, 
and to quantify the global properties of this mode, i.e. excitation energies and 
B(E1) strength, which may also be interesting for the future experimental 
studies. Furthermore, we employ the CRPA model to investigate the role of 
the coupling to the continuum for the proton PDR. We choose to analyse a 
doubly magic nucleus ^^Ni, which is the most proton-rich isotope that has 
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been experimentally discovered [23]. For comparison, we also present results 
for ^^Ni, which lies close to the valley of stability. Being doubly magic, these 
two isotopes can be studied within the CRPA and RRPA models, which do 
not include pairing correlations. Ground state properties of proton-rich nuclei 
around ^^Ni have extensively been studied within the shell-model [24] , Hartree- 
Fock Bogoliubov [25], and Relativistic Hartree-Bogoliubov (RHB) theory [26]. 
Proton drip-line nuclei are characterized by a reduction of the spin-orbit term 
of the effective interaction, outer orbits appear to be very weakly bound, 
and the Fermi energy level may become positive. Due to the presence of the 
Coulomb barrier, loosely bound orbits are stabilised, in contrast to the neutron 
drip-line nuclei where the weakly bound neutron orbits arc more spatially ex- 
tended. The interplay of all these effects in the nuclear ground state will shape 
the properties of the corresponding dipole excitation response. 

In the following, we briefly present the basic theoretical background of the 
Dirac-Hartrcc+RRPA and Skyrme-Hartree-Fock-|-CRPA models and their ef- 
fective interactions. 

In recent years, the relativistic mean-field theory and linear response based 
on density-dependent interactions, turned out to be very successful in studies 
of nuclear ground-state properties and excitation phenomena with a mini- 
mal set of parameters in a fully microscopic way [27]. Within this framework, 
the nucleus is described as a system of Dirac nucleons which interact in a 
relativistic covariant manner by exchange of effective mesons. The model is 
formulated with the Lagrangian density which explicitely includes the density 
dependence in a, a;, and p meson-miclcon vertices. In the present study, we 
employ the density-dependent effective interactions which arc constrained by 
the properties of finite nuclei and nuclear matter: DD-MEl [28], and the new 
interaction DD-ME2 which provides an improved description of the isovec- 
tor dipole response [29]. In the small-amplitude limit, the RRPA equations 
arc derived from the equation of motion for the nucleon density [30]. Both 
in RRPA and CRPA models, we use the same form of the effective isovector 
dipole operator as in Ref. [22]. The RRPA configuration space is constructed 
from particle-hole {ph) pairs composed of the particle states above the Fermi 
level, and hole states in the Fermi sea. In addition, in the relativistic case, one 
also needs to include transitions to unoccupied states from the Dirac sea [27]. 
The resulting RRPA discrete spectra are averaged with the Lorentzian distri- 
bution which includes an arbitrary choice for the width, Tiijipa=^ MeV. The 
Dirac-Hartree+RRPA model is fully self-consistent, i.e. both the equations 
of the ground state, and the residual RPA interaction are derived from the 
same effective Lagrangian. This is an essential property for an accurate de- 
coupling of the spurious center-of-mass motion without need for including any 
additional free parameters. For the present study we use the Dirac-Hartree 
model formulated in the harmonic oscillator basis. Within this approach, the 
particle continuum is represented by a set of discrete states, which are used 
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to construct the RRPA configuration space. In Ref. [26] it has been verified 
that for nuclei towards the proton drip-hne, an expansion in a large oscillator 
basis (N=20) provides sufficiently accurate solutions, in complete agreement 
with the model formulated in the coordinate space. 

The second theoretical framework for the present study is the Skyrmc-Hartree- 
Fock (SHF) plus Contimium-RPA (CRPA) model. The HP equations describ- 
ing the ground state are derived variationally from the Skyrme energy func- 
tional. The ph residual interaction is derived from the same energy functional. 
In the present study, the Coulomb interaction, as well as spin-dependent terms 
are omitted from the residual interaction. The CRPA is formulated in the co- 
ordinate space, and the particle continuum is fully taken into account. The 
transition strength distribution R{E) is continuous by construction. A small 
but finite value of Im£' = r/2 entering the evaluation of the ph Green function 
ensures that bound transitions acquire a finite width and thereby contribute to 
the distribution. More details on the CRPA model can be found in Refs. [31,32] 
and references therein. 

For the purposes of the present study we implement various parameterisations 
of the Skyrme interaction, corresponding to different nuclear-matter proper- 
ties. In particular, they have different (isoscalar) effective mass m*/m and 
isovector effective mass m*^/m. These quantities have influence on the evalu- 
ated properties of the isovector giant dipole resonance (IVGDR) [19,33] and 
the low-energy dipole transitions. In general, interactions with a high effective 
mass are not able to reproduce the IVGDR properties. However, a high m* jm 
may be more appropriate to describe correctly the density of states lying 
close to the Fermi energy, which are relevant for the present study. There- 
fore, we use several different parameterisations of the Skyrme interaction, in 
order to ensure the general validity of our conclusions, at least on a qual- 
itative level. We employ parameterisation MSk7 (m*/m — ml/m — 1.05), 
based on a Hartree-Fock-BCS model [34], and two interactions from the re- 
cent BSk scries, namely BSk8 [35] {m*/m = 0.80, ml/m = 0.87) and BSk2 
{m*/m = 1.04, ml/m = 0.86) [36]. Both BSk interactions were parameterised 
by fitting the values of nuclear masses calculated within the Hartree-Fock- 
Bogohubov method to essentially all the measured ones. We also use the tra- 
ditional parameterisation SkM* {m*/m = 0.786, m^/m = 0.875) [37], which 
was extensively used in previous studies of giant resonances and response of 
exotic nuclei, e.g. Refs. [31,32,38,39]. 

Next, we present the results obtained with the two models. In Fig. 1 we plot 
the ground-state proton and neutron density distributions for ^^Ni and ^^Ni, 
calculated with Dirac-Hartree and SHF models, based on the DD-MEl and 
BSk8 effective interactions respectively. For ^^Ni [Z = 28, N = 28) the pro- 
ton and neutron density distributions are similar in the nuclear interior and 
beyond the surface region the differences completely vanish. However, in the 
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case of ^^Ni {Z = 28, = 20), due to the excess of loosely-bound protons, 
the proton density distribution is considerably extended beyond the neutron 
density distribution. This effect is especially pronounced at radial distances 
r > 2 fm, and it resembles in structure a proton skin. Due to the presence 
of the Coulomb barrier, which tends to localise protons within the nuclear 
interior, the proton skin is not so pronounced an effect as the neutron skin, 
ffowever, some evidence for increasing of the proton-skin thickness in nuclei 
towards the proton drip-line is provided both by theoretical and by experi- 
mental studies [25,40]. 

In Fig. 2 we present results obtained within the fully self-consistent RRPA 
model based on Dirac-Hartree ground state with DD-MEl effective interac- 
tion. In the left panel we display the isovector dipole strength distributions for 
^^Ni and ^^Ni. In the region of the isovector giant dipole resonance (IVGDR), 
the difference between the two distributions is very small and mainly consists 
of small fluctuations of the IVGDR tail. In agreement with the mass depen- 
dence of the giant resonance, the distribution of "^^Ni is only slightly pushed 
to higher energies from the one for ^^Ni. However, in the low-energy region, 
the El strength distributions are rather different: whereas there are no low- 
lying states for ^^Ni, proton-rich ^^Ni is characterized by the appearance of 
a pronounced amount of low-lying transition strength. In order to clarify the 
origin of this strength, in the right panel of Fig. 2 we show the neutron and 
proton transition densities for two characteristic peaks: the low-lying state 
at 7.72 MeV, and the giant resonance state at 18.71 MeV. The transition 
densities of the latter display the dynamics of IVGDR: collective oscillation 
of neutrons against protons. For the low-energy peak, however, the proton 
and neutron transition densities are in phase in the nuclear interior, whereas 
beyond the surface region there are no contributions from the neutrons and 
proton transition density dominates. This type of nuclear dynamics is charac- 
teristic for the proton FDR, where loosely bound protons oscillate against the 
rest of the nucleons [22] , in analogy to the neutron PDR in neutron-rich nuclei. 
The RRPA amplitudes for the El state at 7.72 MeV reveal the structure of 
the proton PDR in detail: the main contribution to the strength of the peak 
consists of transitions from the proton lf7/2 state, which is located at 0.11 
MeV and weakly bound partly due to the presence of the Coulomb barrier. 
The contributions from other transitions are at least an order of magnitude 
smaller. Therefore, the appearance of the low-lying proton PDR strength is 
directly related to the pronounced proton density distributions from Fig. 1. 
Protons from the same loosely-bound orbit contribute to the exotic nuclear 
structure of the ground state, and to the excitation phenomena of the proton 
PDR. The collectivity of the proton PDR peak considerably increases in open 
shell-nuclei, due to the increased number of two-quasiparticle configurations 
composed from many states around the Fermi surface which are, in that case, 
partially occupied [22] . 
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The properties of the low-lying dipole transition strength are strongly sensitive 
on the proton excess. In comparison of the cases of ^^Fe (from Ref. [22]) and 
"^^Ni (DD-MEl interaction is employed in both cases), one can see that the 
GDR peak energy only weakly changes (0.2 MeV) with addition of two more 
protons. The peak energies of the proton PDR mode, however, lowers from 9.4 
MeV towards 7.7 MeV for ^^Fc and ^^Ni, respectively. Obviously, the properties 
of the low-lying dipole transition strength are more sensitive to the variations 
of the nucleon excess than GDR. This type of behaviour indicates the nature of 
the proton PDR mode: as the number of protons increases, the oscillations of 
loosely-bound protons acquire lower frequencies due to their weaker binding. 

By using the Skyrme Hartree-Fock+CRPA model with BSk8 effective interac- 
tion, we repeat the same study of the isovector dipole transition strength for 
^^Ni and ^^Ni isotopes. In the left panel of Fig. 3, we notice that for the two 
Ni isotopes the strength distributions are not very different for E>10 MeV. In 
agreement with the RRPA results, the CRPA low-lying transition strength for 
Ni is strongly enhanced in comparison with the ^^Ni case. In the right panel 
of Fig. 3 we plot the transition densities for a low energy peak at 9.72 MeV 
and the IVGDR at 20.28 MeV. The displayed CRPA transition densities are 
in full agreement with the RRPA calculations, i.e. the high-energy peak corre- 
sponds to the collective IVGDR where protons oscillate versus neutrons, while 
the low-lying transitions reveal the nature of the proton PDR. 

Within the CRPA model, we do not obtain only one pronounced low-energy 
PDR peak for ^^Ni, but rather a smooth continuum, slightly structured around 
9-10 MeV. One of the structures is the 9.72 MeV peak for which the transition 
density was plotted in Fig. 3. This continuum is not an artifact of the small 
smearing parameter used, F = 0.05 MeV. The particle threshold energy is 
E'th = 3.3 MeV. We have evaluated the proton- and neutron-transition densi- 
ties corresponding to various values of excitation energy below £■ = 10 MeV, 
and verified that tlicy are characterized by transitions of a similar nature as the 
9.72 MeV peak. The above discussion remains valid qualitatively when other 
Skyrme interactions are used. More numerical results are presented below. 

In Table 1 we compare the global properties of the proton PDR for ^^Ni, 
obtained using the RRPA model with DD-MEl and DD-ME2 effective inter- 
actions, and the CRPA model with Skyrme interactions BSkS, BSk2, MSk7 
and SkM*. We have calculated the summed low-lying strength mo, the energy- 
weighted strength mi and the centroid energy mi/rriQ for excitation energies 
below 10 MeV, where corresponds to the A;— th moment of the strength 
distribution. Since the CRPA strength distribution is continuous in this re- 
gion, the choice of the cut-off value of PDR region at 10 MeV is somewhat 
arbitrary and the definition of the centroid energy becomes ambiguous. For 
this reason, the respective results are placed inside brackets. In addition, we 
list the relative amount of the low-lying strength mi with respect to the clas- 
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sical Thomas-Reiche-Kuhn sum rule TRK= 14.9 • {NZ/A)e^fm'^ McV. The 
centroid energies of the proton PDR for ^^Ni are obtained around 8 MeV. The 
two different models are in fair agreement for various interactions, exhausting 
for the proton PDR from ^ 0.9% (BSk8) towards 1.6% (DD-MEl, MSk7) 
of the classical TRK sum rule. 

We found that that BSkS and SkM* gave the lowest values of low-lying PDR 
strength. In general, the Skyrme-type interactions with the lower value of 
m*/m, result in weaker low-energy transition strength. The amount of the 
low-lying strength is directly related to the energy e/ of the least bound proton 
single-particle state I/7/2 and the particle threshold energy E^^. For ^^Ni, these 
quantities are as follows: e/ = -0.88 McV and Eth = 4.2 MeV for SkM*, e/ = 
-0.03 MeV and ^Jth = 3.3 MeV for BSk8, e/ = 0.26 MeV and Eth = 3.0 MeV 
for BSk2, and e/ = 0.35 MeV and ^th = 3.0 MeV for MSk7. As the energy 
of 1/7/2 proton state changes to higher values, the corresponding low- lying El 
strength is more enhanced. It is surprising that the Skyrme interaction whose 
effective-mass properties differ the most from the ones of the rclativistic forces, 
namely MSk7, gives the best agreement with the RRPA results. The opposite 
holds for BSkS and SkM*, whose properties differ the least from the ones of the 
relativistic forces. We notice, moreover, that the BSk2 corresponds to almost 
as large an effective mass m* /m as the MSk7, and practically the same i?th, 
yet its predictions for the PDR are closer to the BSk8. Its isovector effective 
mass, though, is low and almost the same as BSk8. Within the present CRPA 
model, a Skyrme interaction cannot provide as much low-lying strength as 
the RRPA model, unless a high value of isovector effective mass is used. It 
is an interesting trend, but it cannot be established only by the study of a 
single nucleus. Finally, there does not seem to exist a correlation between 
the degrees of quantitative agreement of the RRPA and CRPA models on 
the properties of the PDR on one hand, and the e/ value on the other. The 
absolute difference between the value of e/ corresponding to the Skyrme force 
MSk7 (best agreement with RRPA on PDR strength) and the RRPA result, 
0.11 MeV, is almost as large as in the case of the SkM* (the worst agreement). 

In the last row of Table 1, we show an estimate of the PDR width, given by the 



mean deviation of the strength distribution up to 10 MeV, o — J {m^jmi^ — {mi/rrioy. 



Different Skyrme interactions result in similar widths, cr ~ 1.4 MeV. On the 

other side, the width evaluated from the RRPA discrete strength distribution 
with the DD-MEl interaction with a smaller effective mass, results in a = 
0.77 MeV. In the RRPA case it provides only a measure of the fragmentation 
of the low-lying strength, and therefore, it is smaller in comparison with the 
CRPA width which includes the contributions of the continuum. 

In conclusion, we have studied the low-lying dipole response of a represen- 
tative case of a proton drip-line nucleus, namely the doubly magic "^^Ni. We 
have employed two different microscopic models, Dirac-Hartree-|-RRPA, and 
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Skyrme-Hartree-Fock+CRPA. with various effective interactions. Within the 
latter approach, a proper treatment of the particle continuum is included, en- 
abling us to study the relevance of the continuum for the low-lying strength. 
The comparison of El strength distributions for ^^Ni and ^^Ni, and tran- 
sition densities, show that the low-lying dipole strength is a fundamental 
property of the proton-rich nuclei, and it corresponds to the proton pygmy 
dipole resonance, where loosely bound protons vibrate against the approx- 
imately isospin-saturated proton-neutron core. The coupling to the particle 
continuum results in an enhanced width of the proton PDR mode, estimated 
around 1.4 MeV. However, it remains unresolved why the best agreement on 
the proton PDR properties, between the relativistic and nonrelativistic mod- 
els, is obtained for DD-MEl and MSk7 interactions which have quite different 
properties. Whereas DD-MEl represents an advanced density-dependent in- 
teraction appropriate for the studies of both stable and exotic nuclei, MSk7 
has large isoscalar and isovector effective mass, and its properties were recently 
improved in BSk series. A consistent comparison would necessitate a proper 
inclusion of the continuum in the RRPA model, and on the other side, a fully 
self-consistent CRPA model without neglecting of terms in the residual inter- 
action. Nevertheless, we have shown that the two different models agree fairly 
well on the global properties of the proton PDR, and provide a clear theoret- 
ical picture for its underlying nature. We hope that the future experimental 
studies towards the proton drip-line will provide evidence for this exotic mode. 
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Fig. 1. Proton and neutron ground-state density distributions for ^^Ni and ^^Ni 
calculated with the Dirac-Hartree (DD-MEl interaction) and Skyrme-Hartree-Fock 
(BSkS interaction) models. 
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a[MeV] 


0.77 


0.79 


[1.42] 
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[1.48] 


[1.40] 



Table 1 

The low- lying (E<10 MeV) dipole transition strength (mo), energy- weighted low- 
lying strength (mi), its fraction of the classical TRK sum riilc. ccntroid energies 
and widths (cr) for "^^Ni, calculated with RRPA (DD-MEl and DD-ME2 effective 
interactions) and CRPA (BSK8, BSK2, MSk7 and SkM* Skyrme effective interac- 
tions) models. The definition of the centroids and widths in the CRPA model is 
somewhat ambiguous (see text) and therefore the corresponding values are given in 
square brackets. 
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Fig. 2. The Dirac-Hartree-|-RRPA isovector dipole strength distribution for ^^Ni and 
^^Ni (left panel). The proton and neutron transition densities are displayed for the 
low-lying state (E=7.72 MeV) and isovector giant resonance (E=18.71 MeV). The 
DD-MEl effective interaction is employed. The shaded area denotes the difference 
between the low-energy part of the strength distributions of ^^Ni and ^^Ni. 
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Fig. 3. Same as in Fig. 2, but calculated with the Skyrme Hartree-Fock+CRPA 
model by using BSk8 effective interaction. Transition densities are shown for the 
states with E=9.72 and 20.28 MeV. 
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